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ABSTRACT 

We carry out a numerical simulation depicting the effects of a Coronal Mass Ejection (CME) event on 
a close-in giant planet in an extrasolar system. We drive the CME in a similar manner as in simulations 
of space weather events on Earth. The simulation includes the planetary orbital motion, which leads 
to the forming of a comet-like planetary magnetotail which is oriented almost perpendicular to the 
direction of propagation of the CME. The combination of this feature and the fact that the CME does 
not expand much by the time it reaches the planet leads to a unique CME-magnetosphere interaction, 
where the CME itself is highly affected by the presence of the planetary magnetosphere. This change 
in the CME properties throughout the event cannot be estimated by simple, analytic calculations. 
We find that the planet is well-shielded from CME penetration, even for a relatively weak intrinsic 
magnetic field. The planetary angular momentum loss associated with such an event is negligible 
compared to the total planetary angular momentum. We also find that the energy which is deposited 
in the magnetosphere is much higher than in the case of the Earth, and our simulation suggests there 
is a large-scale change in the orientation of the magnetosphere-ionosphere current system during the 
CME event. 

Subject headings: stars: coronae - planet-star interactions - Sun: coronal mass ejections (CMEs) 



1. INTRODUCTION 

Coronal Mass Ejections (CMEs) are sudden releases of 
ionized gas from the stellar surface into interplanetary 
space. These explosions are triggered by a mechanism 
that is not yet fully understood, but it is thought to 
involve slow storage of magnetic energy that is quickly 
released by means of magnetic reconnection. In the so- 
lar case, following the eruption, CMEs are accelerated 
in the corona, some to high terminal speeds of up to 
3000 km s~^, where their magnetic integrity (commonly 
referred to as the "magnetic cloud") is maintained as 
they are propagating to large interplanetary distances 
(for a general review on t he initiation, evolu tion, and 
propagation of CMEs, see Forbes et al. 2006). By the 
time it reaches 1 AU, the interplanetary CME would 
have expanded and reached a size much bigger tha n the 



magnetosphere of the Earth (Forsyth et al. 2006 1. As 



a result, the ambient solar wind conditions m the vicin 
ity of the Earth are replaced by the plasma conditions 
carried by the CME. In some cases, when the orienta- 
tion of the CME magnetic field is opposite to that of 
the Earth, a geomagnetic storm is driven by the energy 
exchange between the magnetic fields. This involves par- 
ticle acceleration and an increase in radiation, as well as 
other phenomena such as increases in ring curren ts, and 



son & Russell 


1995 


phenomena, t 


leir ir 



fort to predict them, have been collect ively known by 
the umbrella ph rase "Space Weather" (Schwenn 2006 
Pulkkinen|[2007| ). 
A large number of exoplanetary systems have been 
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Schneider"1995'; Mayor et al. 2007). Until recently, most 



observed since the mi d-1990's jM avor fc Queloz||1995 



observed cxoplanets were of the "hot-jupiter" type. This 
class of planets are Jupiter-size gas giants, found in close- 
in orbits at distances as small as 0.01 AU from their host 



star (Schneider 1995 Mayor et al. 2003). Planets with 
sizes of Saturn, Uranus, iNeptune, Super-earth, and even 
Earth have been recently observed by new instruments 
that are p articularly designed for plane t searches, such 
as Kepler ( [Borucki et al. 112010 ), CoRo T ( Auvergne et al. 
20091 , and iWCy 7 '([Walker c tal.|2003[ ). Due to the prox- 



imity of hot-jupiters to their host star, it is possible that 
they can be located within the Alfven point, at which 
the stellar wind becomes super- Alfvenic. Therefore, if 
a hot-jupiter has a substantial internal magnetic field, 
it can affect the stellar corona and the star via Star- 
Planet (magnetic) Interaction (SPI). The conse quences 



of SPI have been discussed based on observations ( Shkol 
nik et al.||2003l l 2005a|b[[2008 '; Kashyap ct al."2008 ||tjaar 



et al||2008[ |Pillitter i etaT 2010; Poppcnhacger eFaT 



2010 1 , tEeoretical models (Cunt z eFaTpOOO 



tions 



Lanzal (20081, '2009, 
Prcussc et al. 2006 



Ip et al 

2010), and numerical simula- 
Lipatov ct aL|2005 Johansson 



et al.p009f |Cohen et al.||2009C( ,2010^ . I'he dominant 



scenario tor SPI is a disruption ot the ambient coronal 
magnetic structure by the planetary magnetic field due 
to the planetary orbital motion and via magnetic recon- 
nection, as well as particle acceleration along field lines 
that connect the planet and the star towards the stellar 
surface. 

The proximity to the host star also causes close-in plan- 
ets to suffer from high doses of X-ray and EUV radia- 



mospheric escape rate (e.g. 


Lammer et al. 2003 


Penz 


et al. 


2008 


Yelle et al. 


2008 


Murray-Clay et al. 


2m\. 



In addition, transient CMEs might contribute to 



plane- 



tary erosion, especially if the planetary magnetic field 



2 



is not strong enough to oppose the dynamic pressure 
of the CME. It is commonly assumed that since stel- 
lar flare rat es increase with stellar magnetic activity (e.g. 
Telleschi et al. 2005) and CMEs are associated with flares 



( |Uudel pUU7p 7 the CME rate should be high for active 
stars. Planetary erosion from such host stellar activity 
could be particularly enhanced. 

This issue is particularly important for planets orbit- 
ing M-dwarfs. Since these stars have low luminosities, 
their habitable zone, which covers the planetary dis- 
tances where the surface temperature of a planet allows 



liquid water to exist ( Kasting et al.|199 3|, is located very 
close to the star. Current planet detection methods are 
well-suite d to find such plan ets which can potentially sus- 



tain life (Scalo et al. 2007). The chance of habitability 



might be reduced, however, due to the high l evel of ac- 



tivity of M -dwarfs. [Khodachenko et aL (2007) and Lam 
mer et al. ( ,2007) have studied the ettect ot CME erosion, 
CME induced ion pickup, and extreme EUV radiation 
on the atmospheres of terrestrial planets orbiting an M 
class star. They have used averaged scaled values of stel- 
lar CMEs based on solar CME statistics, as well as other 
scaling laws to characterize the M-type star and the stel- 
lar CME rate. They also determined the magnetospheric 
standoff distance and quantified the planetary magnetic 
moment necessary to shield the planet from CME ero- 
sion. 



The work presente d in Khodachenko et al. ( 2007 ) and 
Lammer et al. ( 2007 ) is based on scaling an Earth- like in- 
teraction between planets and CMEs to a close-in orbit 
scenario. Here we present a full, three-dimensional, time- 
dependent, Magnetohydrodynamic (MHD) simulation of 
the interaction of a CME with the magnetosphere of a 
close-in hot-jupiter. The ultimate goal of our simulation 
is to study how different this interaction is from the case 
of the Earth due to the different magnetospheric struc- 
ture, the interaction taking place partly in a sub-Alfvenic 
region, and the CME not being fully developed by the 
time it reaches the planet. We also use the simulation 
to constrain the CME penetration depth for a particu- 
lar planetary magnetic field strength and investigate the 
loss of planetary angular momentum as a result of the 
interaction. 

We describe our numerical approach, the specifications 
of the planetary magnetosphere, and the CME initiation 
in f|2j We present the simulation results in j|3|and discuss 
their consequences on close-in planets in ^4f^ Finally, we 
summarize our conclusions in S|5] 

2. NUMERICAL SIMULATION 

2.1. Numerical Approach 

In the numerical simulation presented here, we study 
the response of the magnetic field of a close-in planet 
to a CME event using a single-fluid MHD model. The 
simulation includes both the stellar wind, the CME, and 
the planet. The planet is described by boundary values 
for the density and temperature, and by its magnetic 
field. In this work, we do not include a full description 
of the planetary atmosphere and magnetosphere, their 
response to the extreme stellar radiation, or detailed 
magnetospheric/ionospheric electrodynamics. Here we 
present a first step towards a more physical description 
of magnetosphere-CME interaction of close-in planets, 



which could be achieved in principle, by coupling the 
model for the stellar corona with a model for the plane- 
tary m agnetosphere using a numerical framework ( Toth 



et al.|20 05). 

Uur work here comprises three ingredients: 1) the am- 
bient stellar wind; 2) the planet; and 3) the CME. We use 
the BATS-R-US glo bal MHD model and its adaptation 
for the Solar Coro na ( [Powell et al |1999[lT6th et al.|2005"' 



Cohen et al.|2007 ) to perform a simulation that includes 



aU three components. This model has been adapted to 
simulations performe d for stellar coronae. as well a s ex- 
oplanetary systems (Cohen et al. 2009b|c 2010b). In 



addition, this model has been often used to simulate Sun- 



Manchester et al.l|2004, 2008 Lugaz 


et al. 20071 Cohen 


et al. 2008), as well as studying the ev 


olution of CMEs in 


the solar corona (e.g. Liu et al.||2008 


Jacobs et al.||2009 


Cohen et al. 


20U9a 20iUa|). 




'I'he mode 


1 solves the single-fluid set ot conservation 



laws for mass, momentum, magnetic induction, and en- 
ergy (ideal MHD equations): 



dp 

dt 



V • (pu) - 0, 
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dt 
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dt 



V • (uB - Bu) = 0, 



(1) 



pu 



2mo 



(B-B)u-B(B-u) 

^"0 



Ej = pig ■ u). 



in the inertial frame or in the rotating frame (where the 
centrifugal force is added). Here p, u, B, and p are the 
plasma mass density, velocity, magnetic field, and ther- 
mal pressure, respectively, po is the vacuum permeabil- 
ity, g is the gravitational acceleration, and T = 1.5 is 
the ratio of specific heats (polytropic index) at large dis- 
tances from the Sun. Ej is an energy source term that 
depends on the local value of polytropic index. This term 
is responsible for the acceleration of the ambient stellar 
wind, and it is described in detail in Section [2 . 2 1 below . 

2.2. Ambient Stellar Wind 

The ambient stellar wind is obtained in a semi- 
empirical manner based on an observed relation between 
the terminal solar wind speed, Ugw^ and the magnetic 
fiux tube e xpansion factor, /s , from which the wind is 
originated (I Cohen et al. 20071. Once an initial, poten- 
tial (non-MHD) magnetic field distribution is specified 
from a surface magnetic field map, the value of the ex- 
pansion factor for each magnetic field line is calculated, 
hence Usw from each field line can be determined, as well 
as the total kinetic energy per unit mass of the plasma, 
u^^/2. Using Bernoulli integral, this energy can be used 
to specify the stellar surface value of the polytropic in- 
dex, 7o, assuming the boundary values for the density 
and pressure are known: 



1 



Rr?) 



70 Po 
7o - 1 Po ' 



(2) 



with G, A/q, and Rq being the gravitational constant, 
the stellar mass, and stellar radius, respectively. We then 



3 



can define the energy source term: 



P 



P 



1 ■ 



1 r - 1' 



(3) 



where 1 < 7 < F is a radial function of 79. It can be 
seen that for smaller values of 7 (smaller expansion of 
the magnetic flux tube) , the energy source is larger, and 
so is the acceleration of the wind. At large distances from 



the Sun (~ lOi?©), 7 ^ T, and 



0. 



The main advantage of this approach over the use of an 
ad-hoc wind solutions imposed on the magnetic field (the 
Parker wind solution for example; Parker||1958p is that 
here the topology of the win d depends on the ob served 
magnetic field structure (e.g., Phillips et al.||1995 ). This 
model takes into account the dependency of the wind 
structure on the large-scale magnetic topology in the so- 
lar/stellar corona and produces a bi-modal wind with 
regions of slow, denser wind, (which originates from the 
open/closed field boundary where the expansion factor is 
large), and regions of fast, less dense wind (which origi- 
nates fro m open field regions wh ere the expansion factor 
is smafl) ( [McC omas et al.| |2007[ ). 

In order to avoid complexity and isolate the planet- 
CME interaction, we set the stellar magnetic field to be 
a dipole with an equatorial field strength of 2.5G aligned 
with the rotation axis of the star. Other stellar parame- 
ters used here are matched to the observed stellar param- 
eters of HD 189733, with stellar radius, i?+ = O.76i?0, 



stellar mass 

Vt^ = 11.95d (iSchneider 1995 



0.2Mf7), and stellar rotation 



Mayor et al. 2003 



eriod, 

■ The 

ambient stellar wind is solar-like with terminal sp eeds 
ranging between 265 km s~ 
hen et al.|[2009cl |2010b[ ). 



and 800 km s (see Co 



2.3. Modeling the Planet 

In order to model the planet, we impose an additional 
boundary condition in the simulation domain which is 
constrained by the planetary surface density, temp era- 
ture, and magneti c fiel d in a similar manner as in Co- 
[hen etjil.| (|2009c|) and [Cohen et al] (|2010b|). At eacH 
time-step, the coordinates ot this boundary (or plane- 
tary body) are updated based on the planetary orbit. 
Grid cells inside the body are defined as "ghost cells" or 
"boundary cells" and are forced to have the boundary 
values. The solution is updated given the particular set 
of "boundary cells" at a particular time-step. At each 
time step, the planet's location is updated, and bound- 
ary cells that no longer overlap the planet are returned 
to being a "regular cell" and their values are updated ac- 
cording to the MHD solution. This process is illustrated 
in Figure [1] Here, the second boundary condition for the 
planet is set at a semi-major axis of a = S.SRi, and with 
a radius of Rp = 0.2i?^ w 1.5i?,/, where Rj is Jupiter's 
radius. The boundary conditions for the base density 
and temperature are Up = 10*" cm^^ and Tp — 10^ k, 
respectively. 

Here we do not include the planetary gravity since sev- 
eral tests have shown that the effect of the small plan- 
etary mass on the solution is negligible. This is due to 
the fact that the planet is small compared to the simu- 
lation domain, so that the planetary scale height cannot 
be captured as also other physical features, important on 
a planetary scale. The region near the planet is magneti- 
cally dominated, and so stress balance between the mag- 



netosphere and wind/CME does not depend on density 
(through thermal pressure, which requires a more realis- 
tic model for the planet). Therefore, the density gradient 
near the planet is dominated by the difference between 
the boundary density value and the ambient density of 
the stellar corona, as well as the gradient in planetary 
magnetic pressure. This density gradient is more mod- 
erate than that in reality, and it shows, for example, a 
drop of less than order of magnitude between Rp and 
2i?p, while in reality one should expect a drop 2-3 orders 
of magnitude based on the planetary parameters. 

Planetary rotation is omitted; even though it could be 
implemented via the boundary condition for the velocity 
inside the second body (currently, the boundary condi- 
tion for the planetary velocity is of the orbital velocity). 
We choose to ignore rotation in this simulation in order 
not to apply more complexity to a numerical boundary, 
which already has a delicate stability due to its tiny size. 
We assume that close-in planets do not rotate as fast as 
Jupiter due to their sp indown by tidal locking processes 



( Sanchez-Lavega 



2004 



tospheric dynamics is 



and that the planetary magne- 
arth-like (dominated by the solar 
wind) a nd not Jupiter-lik e (dominated by planetary ro- 
tation) ( |Kivelson fc Russen|r995p . 

Since our goal is to quantity the planetary magnetic 
field strength necessary to shield the planet from CME 
events, we study two cases. One with a weak equatorial 
dipole field strength of 0.5 G (Case A hereafter) and one 
with a stronger equatorial dipole field strength of 1 G 
(Case B hereafter). Jupiter's equatorial field is about 
4.3 G. We start the simulation by letting the wind so- 
lution relax to a steady-state with a stationary (tidally 
locked) planet and then turn on the planetary orbital 
motion. The orbital motion is obtained by updating the 
coordinates of the planetary boundary at each time-step. 
This technique and its implications for the dynamics of 
stellar coronae harb oring close-in planets is described in 
Cohen et al. (2011). We allow the simulation to evolve 
tor halt an orbit so that it does not include any pertur- 
bations generated by the initiation of the orbital motion. 
We then use this solution as our initia l co ndition and 
launch the CME as described below (p^. Figure |2] 
shows the distribution of the number density and tem- 
perature at this pre-emption state over the equatorial 
plain of the simulation domain. The relatively high tem- 
perature of ^ 10^ K in the planetary tail (comparing to 
the boundary value of 10^ K) is due to the fact that it 
contains hot coronal plasma that is being trapped inside 
the planetary magnetosphere as it is sweeping through 
the corona. 

2.4. CME Initiation 

We initiate the CME by superimposing an unstable, 
semi-circular flux rope ba sed on the analytical model by 
Titov fc Demoulin ( 1999 ) o n top of the ambient "initial " 
solution described above (j |2.3| (Roussev et al. 2003). 
In simulations of real solar CME events, the flux rope 
properties are matched to fit the observed properties of 
the source active region and its inversion line. The free 
energy provided to the CME is controlled by an addi- 
tional toroidal field in order to produce the observed lin- 
ear speed of the CME. Here, for the sake of definiteness, 
we use the parameters (except for orientation) matched 
to a CME event that occurred on May 13, 2005. This 



4 



was a typical solar Halo CME event which gained a linea r 
speed of about 1700 km (Gopalswamy et al. 2009), 
and is a Solar, Heliospheric, and liNterplanetary hlnvi- 
ronment (shine|^ campaign event. SHINE campaign 
events are picked m order to promote a community effort 
to study a particular CME event from different aspects 
(observations, theory, and numerical modeling). Here we 
investigate the impact of this CME on the exoplanet, and 
in particular, we study the CME penetration depth for a 
particular planetary field strength. Therefore, we choose 
to set the flux-rope orientation to be exactly opposite 
to the planetary dipole orientation. This way, the CME 
penetration is maximized, supported by magnetic recon- 
nection between the CME and the planetary magnetic 
fields. 

2.5. General Simulation Setup 

We run the simulation for 6 hours after the CME on- 
set time. By this time, the CME reaches the edges of 
the simulation domain. Here, we focus on the early part 
of the simulation when the CME reaches the planet, ap- 
proximately 30 minutes after initiation. We omit later 
stages when the the planetary magnetosphere relaxes 
back to its initial state after the CME has passed, and 
in fact, the CME is never turned off. The simulation 
here is done with a model for the stellar corona. There- 
fore, at this point we are unable to study the planetary 
magnetospheric response to the CME in detail. In or- 
der to maintain the CME integrity and preventing nu- 
merical diffusion of the CME magnetic field, we specify 
the non-uniform grid to have very high-resolution around 
the source active region on the stellar surface with Ax = 
10^^ Ri,, and in the space between the star and the planet 
with Ax — 5 • lO^^i?*. The grid size around the planet 
and along its orbit is Ax = 1.1 • ^0~^R^, = 0.055i?p. 

3. RESULTS 
3.1. CME-Planet Interaction 



Figures |3][5] show the number density, temperature, 
and radial speed respectively, on the equatorial plain, 
at t = 00:20, 01:00, 01:20, and 03:00 hours for Case B. 
An animation of these frames, as well as an animation 
of the three-dimensional magnetic field evolution can be 
found in the online version of the paper (Case-B_n.mov, 
CaseB_T.mov, CaseBM.mov, and ESP.CME_3D.mov), 
The large-scale coronal structure looks similar for Case 
A; the differences between the two cases are described 
later as we touch upon the more detailed aspects of the 
results. 

At t = 00:20h, the CME can be seen as it approaches 
the planetary magnetosphere. The CME front is charac- 
terized by a dense, hot plasma. When a CME is propa- 
gating through the solar corona, it expands and acceler- 
ates while pushing the ambient coronal plasma outward 
and interacting with its magnetic field. In our simulation, 
a similar scenario is obtained until the CME reaches the 
nose of the planetary magnetosphere. Then, the CME is 
deflected and adiabatically expands around the magne- 
tosphere front and envelopes it from the sides (noticeable 
at t = 01:00h). The CME front is slowed down and the 
temperature at its front drops. Later on, at t = 01:20h, 

^ SHINE is the NSF branch for solar and heliospheric physics 



the CME reaches the night side of the planetary mag- 
netosphere, while detaching from part of the planetary 
magnetic tail. At this point, the size of the planetary 
magnetosphere, characterized by the low temperature 
bubble near the planet, is significantly decreased. Fi- 
nally, at t = 03:00, the magnetosphere starts to recover 
from the initial impact by the CME, it gets larger, and 
its tail is stretched in the anti-stellar direction. This is 
due to the ambient plasma, which is now dominated by 
the radial propagation and expansion of the CME. 

3.2. CME Penetration Depth 

The cold bubble around the planet is persistent 
throughout the simulations, as seen in Figures |3][5] for 
Case B with a stronger planetary magnetic fieldf indi- 
cating that the CME never reaches the planetary surface 
directly. In order to quantify the penetration depth of 
the CME towards the planetary surface more precisely, 
we extract three spheres around the planet at height of 
h = 0.5, 1, and 2Rp. We calculate the total mass flux, 
Fp — J pUp ■ da, flowing through each surface throughout 
the simulation, where Up is the outward radial velocity in 
the planetary frame of reference, and da is a surface ele- 
ment of the sphere. Since the density at the surface of the 
planet is higher than the density of the ambient coronal 
plasma, the total mass flux over the spheres at the cho- 
sen distances should be positive (a net mass loss from the 
planet). Therefore, a change in the sign of F means that 
the CME has overcome the positive planetary outflow 
and penetrates through the particular sphere, changing 
the net flux through the surface from positive to nega- 
tive. 

In Figure [6) we show the total mass flux through 
the three spheres as a function of time. The fluxes 
are normalized to their value at t = 00:00h, FpQ = 
7 • 10^^^ g cm'^^ s^^. In Case A, the mass flux through 
the sphere at a hight oi h — 2Rp reaches high negative 
values around t — 00:30h. This is a clear evidence for 
CME penetration to that height. However, the negative 
values of the mass flux for h — IRp are about 10 times 
smaller than the values for h — 2Rp, while the mass flux 
measured at h = 0.5Rp remains positive throughout the 
simulation. The mass flux for h — 2Rp increases around 
t = 01:20h, due to the magnetosphere bouncing back 
from the compression by the CME, overshooting its orig- 
inal location and covering this sphere with higher density 
than the original value at that point. The initial mag- 
netospheric compression is the cause for the mass flux 
increase near t = 01:00h for h = 0.5Rp and h = IRp. 
In Case B, the mass fluxes of all three spheres initially 
increase due to the magnetospheric compression, but the 
flux through the h = 2Rp spheres changes its sign around 
t — 01:10h. Nevertheless, the absolute value of the mass 
flux is much lower than in Case A, indicating that while 
the CME can still penetrate to this depth, its impact 
is much weaker. In Case B, the sign of the total fluxes 
for h — O.bRp and h = \Rp remains positive, indicating 
that the CME does not penetrate to a height oih = IRp 
above the planetary surface (^ 18 grid cells). 

3.3. Synthetic Observations of the Space Weather 
Event 

When a CME reaches the Earth during a space weather 
event, it is observed and measured by satellites. Two 
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of these satellites WIND jAcuna et al ||1995[ ) and ACE 
( [Stone et al |1998[ ), are designed to continuously monitor 
the solar wind near Earth. A typical signature in these 
observations appears when a CME with a strong nega- 
tive z component of the magnetic field reaches Earth, z 
is in the coordinate system associated with the Earth's 
magnetic dipole, which points in the positive z direction. 

In our simulation, we can extract the model plasma and 
magnetic field conditions at a particular time and loca- 
tion. Figure [7] shows an extraction of the synthetic data 
from the simulation along the trajectories of two imagi- 
nary satellites which are located at the planetary substel- 
lar point, one at a height ol h = IRp and one at a height 
of /i = 2Rp. Until the CME shock arrival around 00:30h, 
the satellite located at /i = IRp sees an ambient magne- 
tospheric plasma, while the satellite located aX h — 2Rp 
sees coronal plasma. The shock is clearly seen as a jump 
in all parameters and a decrease of the component to 
large negative values. The changes are not sharp as seen 
at real data taken near Earth, since here the differences 
between the ambient plasma and the CME plasma are 
not as marked as at 1 AU. The cold, less dense mate- 
rial seen in Figure [7] between 00:50h and 01:30h is not 
due to the CME cavity, but due to the planetary magne- 
tosphere expanding and bouncing back after the initial 
impact by the CME. Between 01:30h and 02:20h, the 
magnetosphere shrinks again, causing another increase 
in all parameters. After 02:20h, the synthetic data show 
a slow decrease in plasma number density and tempera- 
ture as a result of a steady flow of the CME. We do not 
recover the initial coronal parameters within the time 
domain of the simulation due to the fact that here we 
never turn the CME "off" as occurs naturally in the real 
event. In this late, steady stage of evolution, the size 
of the planetary magnetosphere increases again due to 
the outflow from the mass source at the planetary sur- 
face, and the magnetospheic tail starts to grow as well. 
The dynamic s des cribed above can be seen in the online 



movies (see 



3.11 



4. DISCUSSION 

In our study of a space weather event on a hot Jupiter, 
introduces unique properties of the interaction between 
the planetary magnetosphere of a close-in planet and a 
CME, based on a particular set of parameters. This 
interaction is clearly different than an interaction that 
takes place far from the star, as in the terrestrial case, 
we also discuss the amount of angular momentum the 
planet loses during the event. 

4.1. Properties of a Close-in CME-magnetosphere 
Interaction 

We find the following features of the close-in planet- 
CME interaction. First, the magnetospheric nose is ini- 
tially located not at the subsolar point, but in the di- 



rection of t he orbital motion (also discussed in Cohen 



et al. 2011). Second, the fast planetary orbital motion 
results m a long comet-like magnetosph eric tail ( S chnei-, 
der et al.||1998[ [Schneiter et al.||2007[ ). Therefore, the 
magnetosphere is almost perpendicular in orientation to 
the direction of propagation of the CME, and the CME 
hits the magnetosphere from the side and not on its nose. 
This orientation involves a contact surface between the 



magnetosphere and the CME that is much larger than in 
magnetosphere-CME interaction on Earth. Third, about 
two hours after the CME onset time, the effect of the 
CME on the planetary magnetosphere becomes steady, 
when the magnetospheric orientation is changed from be- 
ing perpendicular to the CME direction of propagation to 
being aligned with it. The magnetospheric nose at this 
stage is roughly at the substellar point, with a steady 
bow shock in front of it. This configuration resembles 
the steady, known magnetospheric configuration in our 
solar system. The initial, perpendicular, configuration 
and the steady, aligned configuration are shown in the 
upper panel of Figure ^ Fourth, by the (short) time 
the CME reaches the planet, it has not had time to be 
fully developed like CMEs in the solar system observed at 
Earth. It does not expand much before it collides with 
the planetary magnetosphere. Hence, it meets a body 
which is comparable in size to it (see bottom-left panel 
of Figure 8]). In the solar system, the scale of CMEs 
is generally much larger than the planetary magneto- 
spheres they meet by the time they propagate to large 
distances in interplanetary space (see bottom-right panel 
of Figure 8|, at least for the terrestrial planets. Even at 
Jupiter, the CME may not be much larger than the mag- 
netosphere, but it can be considered as a plane "wave" 
that impacts the magnetospheric nose without being af- 
fected and modified itself by the interaction. In the case 
of close-in interactions studied here, the CME hits the 
magnetosphere from the side, it is slowed down and cools 
off as it expands around the magnetosphere. The CME 
is in effect split by the planet and ends up surrounding 
the planetary magnetosphere. That the CME is affected 
to such an extent by the magnetosphere is a unique prop- 
erty of the close-in interaction. 

4.2. CME Impact on the Planet and Magnetospheric 
Configuration 

Due t o th e unique interaction properties described 
above (§ 4.1 1, one might expect the impact of the CME 



on the planet to be stronger than in an interaction which 
takes place farther out in interplanetary space. Never- 
theless, our simulation shows that even with a relatively 
weak planetary magnetic field, the CME does not pen- 
etrate more than 0.5Rp above the surface. This means 
that the planetary atmosphere is well-shielded from ero- 
sion by the dynamic pressure of the CME by even a quite 
modest planetary magnetic fields of 0.5 and 1.0 G. With a 
much weaker field, we would expect a more direct hydro- 
dynamic interaction of the CME with thermal pressur e 
of the planetary upper atmosphere ( Lammer et al.|2007 ). 
Such weak field is unlikely though, since the equatorial 
field strength of Jupiter is about 2.5G. 

Here we studied a CME whose magnetic field orien- 
tation is exactly opposite to the planetary field. This 
topology leads to a continuous magnetic reconnection 
between the CME and the planetary fields at the the 
front of the magnetosphere, as well as behind the planet. 
In the case of the Earth, reconnected field lines at the 
dayside are dragged by the CME plasma to the night 
side, where they reconnect again. During the time the 
field lines are dragged from the day to the night side, 
the magnetotail is stretched and energized, and then 
snapped due to the second r econnection. This p rocess 

1995| ). 



drives a geomagnetic storm (Kivelson & Russell 
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In our simulations, the CME completely surrounds the 
magnetosphere between 00:30h and 01:40h, leading to 
a detachment of a significant fraction of the magneto- 
spheric tail. During this disconnection, a blob of plasma 
is snapped towards the planet and is seen to travel ra- 
dially inwards (as shown in Figure^. The average den- 
sity of this plasma blob is 10^ cm" and it is snapped 
through a distance of about O.Si?^, within 10 minutes, so 
the average blob speed is Ubiob ~ 450 km s^^. From 
the simulation, we estimate the volume of this blob to 
be about V « 10^* cm^. _These parameters give a total 
kinetic energy, Ek — ^pU^iob ' ^ ~ 10^^ ergs. This is 
three orders of magnitude higher than the typical energy 
of 10^^ ergs (or 10^^ J) carried by the magnetotail fast 
flows towards the E arth during a geomagneti c substorm 
( Pulkkinen et al.|20 02; Pulkkincn 2007; Tans kanen et al 
2U02| . However, this is about four-hvc orders ot magni- 
tude lower than the EUV energy flux a planet this size 
receives in 24 hours (the time-scae of the event). Nev- 
ertheless, the CME energy is deposited in an impulsive 
manner so that the effect on the planet could be signifi- 
cant due to the geomagnetic effects. 

Figure 10 shows a meridional cut along the planetary 
magnetosphere at t=00:00h, 01:00h, 02:00h, and 03:00h, 
colored with the magnitude of the field-aligned current 
(current parallel to the magnetic field). The initial con- 
figuration resembles the Earth field-aligned current sys- 
tem, with the two tail lobes stretched behind the planet 
separated by a thin current sheet. As the CME hits the 
planet, the two lobes diverge and become more separated 
from ea ch other in a topology that rese mbles "Alfven 
Wings" ( |Drefl et al |1965||Neubauer|1980[ ). At t=03:00h, 
the trailing lobes have inoved back toward each other. 
Alfven wings are Alfven waves generated around a body 
which moves in a magnetized plasma. These waves prop- 
agate away from the body, creating lobes with different 
properties from the surrounding plasma, and where the 
angle of these lobes depends on the Alfvenic Mach num- 
ber. For low Mach numbers, the lobes appear far apart 
while for high Mach numbers, the lobes are close to each 
other behind the body. The Mach number and the topol- 
ogy of the wings/lobes can have a significant effect on the 
magnetospheric response to the CME, the planetary po- 
lar cap potential, and the way ener gy is deposite d from 
the CME into the magnetosphere ( [Ridley | [20071 ). Here 
there is a clear transition of the lobes during the event 
and the effect on the planetary magnetosphere due to 
this may be significant in the above context. However, a 
detailed numerical model for the magnetosphere itself is 
needed in order to study such a magnetospheric response 
to a CME and is beyond the scope of this work. 



4.3. Angular Momentum Loss During the CME Event 

The magnetospheric plasma is confined within the 
magnetosphere, and is moving with the planetary mag- 
netic fields. Therefore, the disconnection of the magne- 
tospheric tail and loss of associated plasma can involve a 
removal of angular momentum from the planet. It is pos- 
sible to compute the total angular momentum associated 
with the planet as the sum of the angular momentum of 
the planet and the angular momentum of the magneto- 
spheric plasma. The angular momentum of the planet 
itself is: 



J It 



r X p 



aMpUp, 



(4) 



where a — 8.8i?* — 4.6-10^'^ cm is the semimajor axis, the 
planetary mass Mp, is assumed to be equal to Jupiter's, 
Mj = 10"^° g, and the planetary angular velocity, Up w 
100 km s~^ (for orbital period of about 2 days at this 
distance) . These parameters give an approximate plane- 
tary angular momentum of about Jp = 5-10^^ g cm^ s~^ . 
Based on the simulation initial result, we estimate the 
total angular momentum of the magnetospheric plasma, 
Jm = J2iPi^i^i X Ui 10^^ g cm^ s~^, where the sum 
is over all the cells in the magnetosphere (excluding the 
planet itself) . These cells are identified by a low temper- 
ature of T < 5 • 10^ K. Later in the simulation, about 
half of the magnetospheric plasma is removed. Taking 
into account that the plasma lost is less dense than that 
which remains attached to the planet, the planet loses 



about Jj\/ = 5 • 10^^ g 



during the event. If there 



were ten such CME events per year, each of which is 
aimed directly at the planet, we obtain an angular mo- 
mentum loss rate of J = 10""'^^ Jp yr~^. This value is 
very small and cannot be significant source of planetary 
angular momentum loss for any plausible CME rate. 

5. SUMMARY 

We carry out an MHD simulation of a space weather 
event on a close-in planet in a manner similar to that 
for the impact of CMEs on Earth. Our simulation 
yields complicated, non-trivial dynamics of a CME- 
magnetosphere interaction, introducing some unique 
properties that do not occur in the solar system. The 
main feature is that the initial orientation of the plane- 
tary magnetosphere, elongated to a comet-like tail by its 
motion in the stellar wind, is almost perpendicular to the 
CME direction of propagation. As a result, the CME is 
modified by the interaction. The resulting change in the 
CME topology and parameters complicates the dynam- 
ics so it is nearly impossible to be described by simple 
analytical relations. This complexity appears even in the 
case of the most simple initial magnetic topology we have 
used here. Despite its proximity to the host star, we find 
that the planet is well shielded from being eroded by 
the CME, even with a relatively weak intrinsic magnetic 
field of 0.5 G. We also find that the planetary angular 
momentum loss associated with a disconnection of part 
of the planetary tail is negligible compared to the total 
planetary angular momentum. Our simulation suggests 
that the planetary magnetosphere can be significantly af- 
fected by the CME event, and that the energization of 
the planetary magnetospheric-ionospheric system might 
be much higher than in the Earth. It also suggests a 
transition in the magnetospheric Alfven wings configura- 
tion during the event, as well as a rotation of the whole 
current system by 90°. However, our simulation cannot 
provide such detailed information about the planetary 
properties; investigation of these aspects of the interac- 
tion requires a detailed numerical model for the planetary 
magnetosphere . 
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Fig. 1. — Top - the change in cells definition as the planet's position is modified illustrated with the actual grid used in the simulation. 
At t=n, cells that are inside the body (colored in yellow) hold the fix boundary values, while cells outside the body (colored in blue) are 
updated via the MHD solution. At t=n+l, the body moves to its new position, so cells that were inside the body are now outside and are 
updated via the MHD solution, while some cells move inside the body and become "body cells" . Bottom - a line extraction of the number 
density across the body at t=n and t=n+l (along the gray line shown on the left). It can be seen that the high boundary value is moving 
with the planetary position. 
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Fig. 2. — The number density (left) and temperature (right) at the initial pre-eruption stage are shown for a slice in the equatorial plane 
in the simulation. The weak planetary field case (Case A) is shown in the top panels and the strong field case (Case B) is shown in the 
bottom panels. The star and the planet are shown as solid spheres. 
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Fig. 5. — Same as Figurc|3]but for radial speed away from the star. 
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Fig. 6. — Total mass flux through spherical surfaces with radii r = 0.5Rp (black curve), r = IRp (blue curve), and r = 2Rp (red curve) 
above the planetary surface, and normalized to the initial value for Case A (top; weak planetary field) and Case B (bottom; strong planetary 
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Fig. 7. — Synthetic stellar wind values extracted from the simulation, mimicing measurements taken by imaginary satellites assumed to 
be located at the subsolar point at a height of IRp (solid line) and 2Rp (dashed line) above the planetary surface. Data are shown for Case 
A (blue) and Case B (red). 



15 




n [cm ] 
6.0E4-05 
3.1E-h05 
1.9E+05 
1.1E*05 
7.1E-h04 
4.3E+04 
2.7E*04 
1.6E-h04 
1.0E+04 — 




T = 00:0000 



X[R] 



T = 06:00:00 



X[R] 




50 
Y[Rs] 



Fig. 8. — Top: number density in the equatorial plane near the planet at t=00:00h (left) and t = 03:00h. Bottom: comparison of the CME 
and the planetary magnetosphere sizes at different distances from the star - left panel shows CME field lines (red), magnetospheric field 
lines (cyan), and an isosurface of Ur = 1500 km s~^, which represents the CME fro nt, and right panel shows , for comparison, a meridional 
cut colored with contours of U from a simulation of a solar CME at 1 AU (taken from |Manchester et al.| | |2004| ). The black ellipse represents 
a tentative size of the Earth's magnetosphere, while the yellow ellipse represents a tentative size ot Jupiter's magnetosphere. 
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Fig. 9. — Number density values, as in the top panels of Figure [To] showing magnetotail plasma that gets 'snapped' towards the planet 
between 01:20h (left) and 01:30h (right). The arrows represent the direction of movement of the plasma. 
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